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Halanaerobium praevalens Zeikus et al. 1984 is the type species of the genus Halanaero- 
bium, which in turn is the type genus of the family Halanaerobiaceae. The species is of inter- 
est because it is able to reduce a variety of nitro-substituted aromatic compounds at a high 
rate, and because of its ability to degrade organic pollutants. The strain is also of interest be- 
cause it functions as a hydrolytic bacterium, fermenting complex organic matter and produc- 
ing intermediary metabolites for other trophic groups such as sulfate-reducing and methano- 
genic bacteria. It is further reported as being involved in carbon removal in the Great Salt 
Lake, its source of isolation. This is the first completed genome sequence of a representative 
of the genus Halanaerobium and the second genome sequence from a type strain of the fami- 
ly Halanaerobiaceae. The 2,309,262 bp long genome with its 2,110 protein-coding and 70 
RNA genes is a part of the Genomic Encyclopedia of Bacteria and Archaea project. 



Introduction 



Strain GSL T (= DSM 2228 =ATCC 33744) is the type 
strain of the species Halanaerobium praevalens, 
which is the type species of its genus Halanaero- 
bium [1]. Originally described as Haloanaerobium 
[2], the name was later changed to Halanaerobium 
to conform with rule 61 of the Bacteriological Code 
[3]. The genus currently consists of nine validly 
named species [4]. The genus name is derived from 
the Latinized Greek word hals; halos meaning salt, 
the Latinized Greek word an- meaning not, the La- 
tinized Greek word aer meaning air and the Lati- 



nized Greek word bios meaning life, yielding the 
Neo-Latin word 'Halanaerobium' meaning 'salt or- 
ganism which grows in the absence of air' [4]. The 
species epithet is derived from the Latin word 
'praevalens' (very powerful, very strong, here pre- 
valent) [4]. Strain GSL T was isolated from the 
hypersaline surface sediments of Great Salt Lake, 
Utah, USA [2]. Further strains of H. praevalens have 
been isolated from canned salted Swedish fer- 
mented herrings referred to as Surstromming [5] 
and probably also from the Red Sea [6]. Other 
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members of the genus have been isolated also from 
high salt environments distributed worldwide [7- 
10]. The enzymatic activities of the fatty acid syn- 
thetase complex and the D-BAPA (N'-benzoyl- 
arginine-p-nitroanilide)-ase of H. praevalens have 
been studied in more detail [11,12]. Here we 
present a summary classification and a set of fea- 
tures for H. praevalens strain GSL T , together with 
the description of the complete genome sequencing 
and annotation. 

Classification and features 

A representative genomic 16S rRNA sequence of H. 
praevalens was compared using NCBI BLAST under 
default values (e.g., considering only the best 250 
hits) with the most recent release of the Greengenes 
database [13] and the relative frequencies, weighted 
by BLAST scores, of taxa and keywords (reduced to 
their stem [14]) were determined. The five most 
frequent genera were Halanaerobium (81.9%), Halo- 
thermothrix (7.8%), Halanaerobacter (2.7%), Aceto- 
halobium (2.3%) and Natroniella (1.9%). Regarding 
hits to sequences from other members of the genus, 
the average identity within HSPs (high-scoring seg- 
ment pairs) was 97.8%, whereas the average cover- 
age by HSPs was 96.3%. The species yielding the 
highest score was Halanaerobium saccharolyticum. 
(Note that the Greengenes database uses the INSDC 
(= EMBL/NCBI/DDBJ) annotation, which is not an 
authoritative source for nomenclature or classifica- 
tion.) The five most frequent keywords within the 
labels of environmental samples which yielded hits 
were 'microbi' (9.4%), 'hypersalin' (9.1%), 'mat' 
(8.6%), 'len, miniprim, new, view, world' (8.5%) and 
'food' (3.4%). The single most frequent keyword 
within the labels of environmental samples which 
yielded hits of a higher score than the highest scor- 
ing species was 'hypersalin, len, mat, microbi, mini- 
prim, new, view, world' (12.5%). These key words 
are in line with the ecology and the niche from 
where strains of H. praevalens have been isolated. 

Figure 1 shows the phylogenetic neighborhood of 
H. praevalens GSL T in a 16S rRNA gene based tree. 
The sequences of the four 16S rRNA gene copies in 
the genome differ from each other by up to five 
nucleotides, and differ by up to five nucleotides 
from the previously published 16S rRNA gene se- 
quence (AB022034). 

The cells of strain GSL T are straight and rod-shaped 
(1 x 2.4 |im) (Figure 2) when grown to the mid-log 
phase at 37 Q C on CS medium containing 12.5% 
NaCl and 0.5% glucose [2]. When grown at higher 



NaCl concentrations (> 20%) the cells appear gra- 
nulated and shorter in length [2]. Single colonies 
were reported as white to translucent in color, 0.5- 
2.0 mm in diameter, and glistened, when grown on 
agar plates containing CS medium, 12.5% NaCl, and 
0.5% glucose [2]. H. praevalens cells stain Gram- 
negative [2] and electron microscopy in thin sec- 
tion revealed architectural features typical of 
Gram-negative bacteria [2]. However, the positive 
D-BAPA-ase reaction [11] confirms its phylogenetic 
affiliation to the endospore-forming firmicutes 
(Table 1). In this respect, H. praevalens is able to 
hydrolyze only the D- but not the L- isomer of N'- 
benzoyl-arginine-p-nitroanilide (BAPA) [11]. The 
activity of D-BAPA-ase was highest at low NaCl 
concentration (100 mM) and completely inhibited 
at NaCl concentration equal or higher than 1.0 M 
(~12%) [11]. Strain GSL T is described to be non- 
motile, although many flagellar genes have been 
identified in the genome (see below). Other isolates 
of H. praevalens were described as motile [5,6], as 
were other members of the genus [7,9,10], suggest- 
ing strain GSL T is atypical with regard to motility. 
The organism is a strictly anaerobic chemoorgano- 
troph [2]. It grows at NaCl concentrations between 
2% and 30%, with optimal growth at approximate- 
ly 13% [2]. The doubling time is 4 h at 12.5% NaCl 
and 7 h at 25% NaCl in complex CS medium [2]. 
The temperature range for growth ranges from 5 Q C 
to 60 Q C, with an optimum at 37 Q C [2]. The pH range 
for growth is between 6.0 and 9.0, with an opti- 
mum at pH between 7.0 and 7.4 [2]. Strain GSL T is 
able to utilize carbohydrates (including pectin and 
iV-acetylglucosamine), amino acids, yeast extract, 
and trypticase; the two latter serving as carbon and 
energy sources on complex medium [2]. The fer- 
mentation of glucose yielded butyrate, acetate, 
propionate, H2, and CO2 as major products [2]. Also, 
fructose, D-mannose and maltose are utilized and 
methionine is transformed to methylmercaptan [2]. 
Penicillin, tetracycline, cycloserine, chlorampheni- 
col (each at 100 ng/ml culture) or sodium azide 
(500 |ig/ml) completely inhibit the growth of H. 
praevalens [2]. Strain GSL T was also able to degrade 
nitro-substituted aromatic compounds such as ni- 
trobenzene, o-nitrophenol, m-nitrophenol, p- 
nitrophenol, 2,4-dinitrophenol, and 2,4-dinitroaniline 
[34]. The fatty acid synthetase of H. praevalens is 
only slightly inhibited at 17.5% and was the first 
reported to be active in the presence of high salt 
concentrations [12]. H. praevalens was reported to 
be involved in carbon sequestration in the Great 
Salt Lake [35], since it is present in the sediments of 
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this lake in high numbers (> 10 8 cells/ml) [2,36]. H. 
praevalens regulates its internal osmotic pressure 
by the accumulation of salts (Na + , K + , CI ) rather 
than by compatible solutes [36]. High concentra- 
tions of these salts were measured inside the cells, 
in sufficient concentration to be isotonic or hyper- 
tonic with the medium [37]. Thiosulfate reduction 



and rhodanese-like enzyme (thiosulfatexyanide 
sulfur-transferase) activities also tested positive in 
strain GSL T [8]. Early in 1987, Matheson et al. [38] 
established the primary structure of the ribosomal 
A-protein of the strain GSL T , which is the equivalent 
to the ribosomal protein L12 from Escherichia coli. 



Halanaerobium praevalens (IMG2503336955) " 
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Figure 1. Phylogenetic tree highlighting the position of H. praevalens GSL T relative to the other type 
strains within the family Halanaerobiaceae. The tree was inferred from 1,460 aligned characters [15,16] 
of the 1 6S rRNA gene sequence under the maximum likelihood criterion [1 7] and rooted in accordance 
with the current taxonomy. The branches are scaled in terms of the expected number of substitutions per 
site. Numbers to the right of bifurcations are support values from 1,000 bootstrap replicates [18] if larger 
than 60%. Lineages with type strain genome sequencing projects registered in GOLD [19] are labeled 
with an asterisk, published genomes with two asterisks [20]. 
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Figure 2. Scanning electron micrograph of H. praevalens GSL T 
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Table 1. Classification and general features of H. praevalens GSL T according to the MIGS recommendations 



[21] and the NamesforLife database [22]. 


1 11 f* C 1 

MIGS ID 


Property 


Term 


Evidence code 






Domain Bacteria 


TAS 


[23] 






Phylum Firmicutes 


TAS 


[24,25] 






Class Clostridia 


TAS 


[26,27] 




Current classification 


Order Halanaerobiales 


TAS 


[28-30] 




Family Halanaerobiaceae 


TAS 


[1,31] 






r~ ill 1 ■ 

Genus Halanaerobium 


TAS 


[1-3,30] 






Species Halanaerobium praevalens 


TAC 

1 AS 


r 1 ti 

[1,2] 






Type strain GSL 


1 AS 


[2] 




Gram stain 


negative 


TAS 


[2] 




Cell shape 


L 1 1. [ 

straight rods 


TAS 


[2] 




Motility 


non-motile 


TAS 


[2] 




Sporulation 


none 


TAS 


[2] 




Temperature range 


above 5 C and below 60 C 


TAS 


[2] 




Optimum temperature 


37°C 


TAS 


[2] 




c 1 ■ ■ i 

Salinity 


2%-30% NaCI, optimum at 13% 


TAS 


[2] 


K A 1 f C 11 

MICjS-22 


Oxygen requirement 


strictly anaerobic 


TAS 


[2] 




Carbon source 


yeast extract, trypticase 


TAS 


[2] 




Energy metabolism 


1 . 1 

chemoorganotroph 


TAS 


[2] 


MIGS-6 


Habitat 


saline environments 


TAS 


[2] 


MIGS-1 5 


Biotic relationship 


not reported 


NAS 




MIGS-14 


Pathogenicity 


not reported 


NAS 






Biosafety level 


1 


TAS 


[32] 




Isolation 


surface sediments of a saline lake 


TAS 


[2] 


MIGS-4 


Geographic location 


Great Salt Lake, Utah, USA 


TAS 


[2] 


MIGS-5 


Sample collection time 


between August 1 979 and August 1 980 


TAS 


[2] 


MIGS-4.1 


Latitude 


41.15 


NAS 




MIGS-4.2 


Longitude 


-112.67 


NAS 




MIGS-4.3 


Depth 


10 m bottom sediment 


TAS 


[2] 


MIGS-4.4 


Altitude 


1 .755 m above sea level 


NAS 





Evidence codes - IDA: Inferred from Direct Assay (first time in publication); TAS: Traceable Author 
Statement (i.e., a direct report exists in the literature); NAS: Non-traceable Author Statement (i.e., not 
directly observed for the living, isolated sample, but based on a generally accepted property for the 
species, or anecdotal evidence). These evidence codes are from of the Gene Ontology project [33]. If 
the evidence code is IDA, the property was directly observed by one of the authors or an expert men- 
tioned in the acknowledgements. 

Chemotaxonomy 

When grown on CS medium, at 5% NaCI, more li- 
pids are produced than at 25% NaCI (3.74% and 
2.54% of the dry weight of the organism, respec- 
tively) [2]. At 5% NaCI, the fractions of glycolipids 
and phospholipids are 46.9% and 44.5% of the 
total lipids, respectively. At 25%, the proportion 
changes in favor of phospholipids (49.1%), whe- 
reas glycolipids decrease (43.0%) [2]. The glycoli- 
pids consist of a single component diacylglycerol 
derivative, while the phospholipids consist mainly 
of cardiolipin (CL), phosphatidyl glycerol (PG), 



and three minor unidentified constituents [2]. 
When grown on CS medium, at 5%, the major fatty 
acids are Ci 4: o (49.3%), Ci 6: i (31.3%) and Ci 6: o 
(11.4%). At 25%, these fractions change to 36.8%, 
39%, and 22.7%, respectively. Similar though 
more detailed results on the fatty acid composi- 
tion have been reported recently [6]. 
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Genome sequencing and annotation 

Genome project history 

This organism was selected for sequencing on the 
basis of its phylogenetic position [39], and is part 
of the Genomic Encyclopedia of Bacteria and Arc- 
haea project [40]. The genome project is depo- 



sited in the Genome On Line Database [19] and the 
complete genome sequence is deposited in Gen- 
Bank. Sequencing, finishing and annotation were 
performed by the DOE Joint Genome Institute 
(JGI). A summary of the project information is 
shown in Table 2. 



Table 2. Genome sequencing project information 



MIGS ID 


Property 


Term 


MIGS-31 


Finishing quality 


Finished 


MIGS-28 


Libraries used 


tI 1 • 1 ' a r- a l 1 11*1 

Three genomic libraries: one 454 pyrosequence standard library, 


one 454 PE library (1 3 kb insert size), one lllumina library 


MIGS-29 


Sequencing platforms 


lllumina GAM, 454 GS FLX Titanium 


MIGS-31. 2 


Sequencing coverage 


201.2 x lllumina; 174.2 x pyrosequence 


MIGS-30 


Assemblers 


Newbler version 2.3, Velvet, phrap 


MIGS-32 


Gene calling method 


Prodigal 1.4, GenePRIMP 




INSDC ID 


CP002175 




Genbank Date of Release 


October 21, 2010 




GOLD ID 


Gc01415 




NCBI project ID 


32591 




Database: IMG-GEBA 


2503283011 


MIGS-13 


Source material identifier 


DSM2228 




Project relevance 


Tree of Life, GEBA 



Growth conditions and DNA isolation 

H. praevalens GSL T , DSM 2228, was grown anae- 
robically in DSMZ medium 210 [' Haloanaerobium' 
medium) [41] at 30-37°C. DNA was isolated from 
0.5-1 g of cell paste using MasterPure Gram- 
positive DNA purification kit (Epicentre 
MGP04100] following the standard protocol as 
recommended by the manufacturer, with modifi- 
cation st/DL for cell lysis as described in Wu et al. 
[40]. DNA is available through the DNA Bank Net- 
work [42]. 

Genome sequencing and assembly 

The genome was sequenced using a combination 
of lllumina and 454 sequencing platforms. All 
general aspects of library construction and se- 
quencing can be found at the JGI website [43]. Py- 
rosequencing reads were assembled using the 
Newbler assembler (Roche). The initial Newbler 
assembly, consisting of 85 contigs in 31 scaffolds, 
was converted into a phrap [44] assembly by mak- 
ing fake reads from the consensus, to collect the 
read pairs in the 454 paired end library. lllumina 
sequencing data (360 Mb) was assembled with 
Velvet [45] and the consensus sequences were 
shredded into 1.5 kb overlapped fake reads and 
assembled together with the 454 data. The 454 



draft assembly was based on 401.6 Mb 454 draft 
data and all of the 454 paired end data. Newbler 
parameters are -consed -a 50 -1 350 -g -m -ml 20. 
The Phred/Phrap/Consed software package [44] 
was used for sequence assembly and quality as- 
sessment in the subsequent finishing process. Af- 
ter the shotgun stage, reads were assembled with 
parallel phrap (High Performance Software, LLC). 
Possible mis-assemblies were corrected with ga- 
pResolution [43], Dupfinisher [46], or sequencing 
cloned bridging PCR fragments with subcloning. 
Gaps between contigs were closed by editing in 
Consed, by PCR and by Bubble PCR primer walks 
(J.-F. Chang, unpublished). A total of 417 addition- 
al reactions and ten shatter libraries were neces- 
sary to close gaps and to raise the quality of the 
finished sequence. lllumina reads were also used 
to correct potential base errors and increase con- 
sensus quality using a software Polisher devel- 
oped at JGI [47]. The error rate of the completed 
genome sequence is less than 1 in 100,000. To- 
gether, the combination of the lllumina and 454 
sequencing platforms provided 375.4 x coverage 
of the genome. The final assembly contained 
838,597 pyrosequence and 12,903,210 lllumina 
reads. 
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Genome annotation 

Genes were identified using Prodigal [48] as part 
of the Oak Ridge National Laboratory genome an- 
notation pipeline, followed by a round of manual 
curation using the JGI GenePRIMP pipeline [49]. 
The predicted CDSs were translated and used to 
search the National Center for Biotechnology In- 
formation (NCBI) non-redundant database, Uni- 
Prot, TIGR-Fam, Pfam, PRIAM, KEGG, COG, and In- 
terPro databases. Additional gene prediction anal- 
ysis and functional annotation was performed 
within the Integrated Microbial Genomes - Expert 
Review (IMG-ER) platform [50]. 



Genome properties 

The genome consists of a 2,309,262 bp long chro- 
mosome with a G+C content of 30.3% (Figure 3 
and Table 3). Of the 2,180 genes predicted, 2,110 
were protein-coding genes, and 70 RNAs; 42 
pseudogenes were also identified. The majority of 
the protein-coding genes (77.7%) were assigned 
with a putative function while the remaining ones 
were annotated as hypothetical proteins. The dis- 
tribution of genes into COGs functional categories 
is presented in Table 4. 



2300001 

2200001 100001 




1200001 1100001 



Figure 3. Graphical circular map of the chromosome. From outside to the center: Genes on forward strand 
(color by COG categories), Genes on reverse strand (color by COG categories), RNA genes (tRNAs green, 
rRNAs red, other RNAs black), GC content, GC skew. 
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Table 3. Genome Statistics 



Attribute 


Value 


% of Total 


Genome size (bp) 


2,309,262 


100.00% 


DNA coding region (bp) 


2,059,925 


89.20% 


DNA G+C content (bp) 


699,559 


30.29% 


Number of replicons 


1 




Fxtr^rhrornnsorna 1 elements 


o 




Total genes 


2,180 


100.00% 


RNA penes 


70 


3.21 % 


rRNA operons 


A 

4 




Protein-coding genes 


z, 1 1 U 


7Q0/ 


rseuuu genes 


HA 


I .yj to 


(^finf^Q \A/itn function nrAniptinn 

v_J V„ 1 1 C3 VV llll IUIILLIWII Ul CUILLIUI 1 


1 694 


77 71 % 

/ / . / I /O 


Genes in paralog clusters 


310 


14.22% 


Genes assigned to COGs 


1,760 


80.73% 


Genes assigned Pfam domains 


1,860 


85.32% 


Genes with signal peptides 


446 


20.46% 


Genes with transmembrane helices 


582 


26.70% 


CRISPR repeats 


1 





Table 4. Number of genes associated with the general COG functional categories 



Code 


value 


%age 


Description 


J 


133 


6.9 


Translation, ribosomal structure and biogenesis 


A 


0 


0.0 


RNA processing and modification 


K 


124 


6.4 


Transcription 


L 


115 


6.0 


Replication, recombination and repair 


B 


1 


0.0 


Chromatin structure and dynamics 


D 


24 


1.2 


Cell cycle control, cell division, chromosome partitioning 


Y 


0 


0.0 


Nuclear structure 


V 


18 


0.9 


Defense mechanisms 


T 


126 


6.5 


Signal transduction mechanisms 


M 


106 


5.5 


Cell wall/membrane/envelope biogenesis 


N 


76 


3.9 


Cell motility 


Z 


0 


0.0 


Cytoskeleton 


W 


0 


0.0 


Extracellular structures 


u 


44 


2.3 


Intracellular trafficking, secretion, and vesicular transport 


o 


66 


3.4 


Posttranslational modification, protein turnover, chaperones 


c 


123 


6.4 


Energy production and conversion 


G 


150 


7.8 


Carbohydrate transport and metabolism 


E 


136 


7.0 


Amino acid transport and metabolism 


F 


74 


3.8 


Nucleotide transport and metabolism 


H 


80 


4.1 


Coenzyme transport and metabolism 


I 


45 


2.3 


Lipid transport and metabolism 


P 


121 


6.3 


Inorganic ion transport and metabolism 


Q 


18 


0.9 


Secondary metabolites biosynthesis, transport and catabolism 


R 


186 


9.6 


General function prediction only 


S 


168 


8.7 


Function unknown 




420 


19.3 


Not in COGs 
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